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ABSTRACT 

We present the first collective evidence that Fermi - detected jets of high kinetic power (Akin) are 
dominated by inverse Compton emission from upscattered external photons. Using a sample with 
a broad range in orientation angle, including radio galaxies and blazars, we find that very high 
power sources (Akin > 10 45 ' 5 erg s _1 ) show a significant increase in the ratio of inverse Compton 
to synchrotron power (Compton dominance) with decreasing orientation angle, as measured by the 
radio core dominance and confirmed by the distribution of superluminal speeds. This increase is 
consistent with beaming expectations for external Compton (EC) emission, but not for synchrotron- 
self Compton (SSC) emission. For the lowest power jets (Akin < 10 43 5 erg s _1 ), no trend between 
Compton and radio core dominance is found, consistent with SSC. Importantly, the EC trend is 
not seen for moderately high power flat spectrum radio quasars with strong external photon fields. 
Coupled with the evidence that jet power is linked to the jet speed (Kharb et al. 2010), this Ending 
suggests that external photon fields become the dominant source of seed photons in the jet comoving 
frame only for the faster and therefore more powerful jets. 

Subject headings: galaxies: active — quasars: general — radiation mechanisms: non-thermal 


1. INTRODUCTION 

The relativistic jets of radio-loud active galactic nu- 
clei (AGN) are copious gamma-ray emitters, as was first 
discovered by EGRET (Hartman et al. 1992), and con- 
firmed by the Fermi large area telescope (LAT), which 
has associated over 800 sources with radio-loud AGN in 
the second catalog (2LAC; Ackermann et al. 2011). Most 
of these are blazars, seen with the jet axis along the line 
of sight, though Fermi has also detected the jets of sev- 
eral radio galaxies (RG; Abdo et al. 2010; Kataoka et al. 
2011), which are misaligned blazars under the standard 
unification scheme (Urry & Padovani 1995). 

The lower-energy peak in the jet spectrum is well- 
understood as synchrotron emission from relativistic 
electrons in the jet. The high-energy component, peak- 
ing from X-ray to TeV energies, is attributed to pho- 
tons upscattered by the same relativistic electrons to 
higher energies via the inverse Compton (IC) pro- 
cess (see Bottcher 2007 for a review). These pho- 
tons could arise from the jet synchrotron emission (syn- 
chrotron self-Compton emission, SSC; Maraschi et al. 
1992; Marscher & Travis 1996) or from external sources 
such as the accretion disk (Dernier et al. 1992), broad- 
line region (BLR; Sikora et al. 1994), or molecular torus 
(MT; Blazejowski et al. 2000; Sikora et al. 2009), i.e. ex- 
ternal Compton (EC) emission. Identifying the IC emis- 
sion mechanism is a diagnostic for the location of the 
gamma-ray emitting region, a currently open issue (e.g. 
Agudo et al. 2011). However, the spectral energy dis- 
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tributions (SEDs) of individual sources are rarely suf- 
ficiently constraining of the IC mechanism due to the 
number of free parameters (Sikora et al. 1997). 

The beaming pattern (how apparent luminosity 
changes with orientation) is different for EC and SSC 
emission. Thus, a collective study using a sample of 
sources at different orientations can be used to identify 
the gamma-ray emission mechanism. In this letter we 
discuss the effect of jet power and orientation on the ob- 
served IC power in view of the recent suggestion of a 
dichotomy in radio- loud AGN (§2). We then show that 
the collective beaming pattern for a sub-set of high-power 
sources appears to support EC models rather than SSC 
for this population (§3), and discuss the implications for 
the gamma-ray emission region. 

2. A DICHOTOMY IN RADIO-LOUD AGN 
2.1. The Synchrotron Plane 

Using a large sample of jets, we recently found ev- 
idence that radio-loud AGN form two populations in 
the plane of synchrotron peak luminosity (A PjSyn ) ver- 
sus peak frequency (^ PlS y n ; Meyer et al. 2011, hereafter 
Mil). A population of ‘weak’ jets consists of sources 
with low A P)Syn which appear to be most aligned at high 
z/ PiSyn (-10 1 ' Hz) and trace out a shallow track on the 
plane as they become less aligned (dropping more in 
^p,syn than A PjSyn ). A separate population of ‘strong’ 
jets with higher A PjSyn and t'p.syn < 10 15 Hz appear to 
drop rapidly in luminosity with misalignment (as mea- 
sured by radio core dominance). Importantly, the weak 
jets exhibit jet kinetic powers below 10 44,5 erg s _1 while 
all sources with Akin > 10 445 erg s _1 are on the strong- 
jet branch. The weak/strong divide in the synchrotron 
plane, which we associate with the morphological di- 
chotomy in Fanaroff- Riley (FR) type 1 and 2 RG, may 
also be mapped to a critical transition in accretion effi- 
ciency (Ghiscllini & Celotti 2001; Ghisellini et al. 2009a; 
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Fig. 1. This plot of total LAT-band luminosity versus gamma- 
ray energy index a (approximating T p ,ic versus ^ p ,ic) is the high- 
energy analog of the synchrotron plane discussed in Mil. Blazars 
have been divided based on the luminosity and location of their 
synchrotron peaks into strong-jet and weak-jet blazars as described 
in the text. Those without ^ P)S yn and L P)S yn are ‘unclassified’. 
FR2 RG appear to have similar z/p, S yn (similar a) to the strong- 
jet blazars while FRI RG have much lower values of zv P)S yn (higher 
a) relative to the weak-jet blazars with which they are associated. 
The FRI would be detectable by Fermi even at much harder a (see 
Ackermann et al. 2011, Figure 15). The region at upper right re- 
mains empty despite recent claims of high-power, high-peak sources 
(shown as black stars, Padovani et al. 2012), matching what has 
been found in the synchrotron plane. 

Georganopoulos et al. 2011). 


2.2. The Blazar Envelope at High Energies 

In light of the dichotomy discussed above, we present in 
Figure 1 the total LAT-band luminosity 5 (a proxy for the 
IC peak luminosity, F Pi ic), versus the LAT-band energy 
index a (a proxy for peak frequency, i-'p.ic) f° r the entire 
Fermi-detected sample of radio-loud AGN of known red- 
shift (data taken from the 2LAC). Sources with higher 
F'p jc will have harder (smaller) values of a, thus in the 
figure r'pjc increases to the right. 

An empty region at upper right appears analogous to 
that seen in the synchrotron plane. Gionnni et al. (2012) 
have suggested that this is a selection effect due to a lack 
of redshifts for high-frequency-peaked, liigh-luminosity 
sources, and Padovani et al. (2012) discuss four such can- 
didates. However, these candidates do not cross into the 
upper-right region in Figure 1 (black stars). These inter- 
esting sources appear to exhibit v PiSyn at or above ~ 10 15 
Hz; however, the SED sampling is sparse, and it is dif- 
ficult to rule out an alternative explanation for the soft 
X-ray spectra such as an extra emission component (see 
e.g., the case of BL Lacertae; Raiteri et al. 2010). Im- 
portantly, for the several hundred Fermi sources lacking 
identifications or redshifts, nearly all (>99%) have soft 

5 tloOMeV— lOOGeV = 47rd^ S-y /(I + z) 1- “ for power-law sources, 
where S-y is the catalog integrated energy flux and a = F — I using 
the published photon index IF For sources fit with Log-Parabolic 
spectra, the integrated luminosity is calculated from the 2LAC 
fitted values with a k-correction. 


spectra (a > 1) or would require redshifts > 2 or higher 
to cross into this region, making it unlikely to be empty 
due to selection effects. 

The blazars in Figure 1 have been divided in an ap- 
proximate way based on their synchrotron spectra into 
strong and weak jets. We take the former to be those 
defined by F P)Syn > 10 45 erg s _1 , v P , sy n < 10 15 Hz, as 
well as flat spectrum radio quasar (FSRQ) type sources 
outside this area. For the weak jets, we take all sources 
outside this area, minus any FSRQ. While very rare cases 
of FRl-like (e.g., weak) FSRQs exist, their occurrence is 
negligible in the bright Fermi-detected population con- 
sidered here (see e.g., Kharb et al. 2010, and references 
therein). The eleven Fermi-detected RG are also shown 
(data from Abdo et al. 2010). The FRI in Figure 1 have 
much lower i/ p jc (larger a) compared with the weak-jet 
sources, while there is little difference between a for FR2 
and strong jets, tentatively matching the different mis- 
alignment paths for these populations that were found in 
the synchrotron plane. 

2.3. The Importance of Jet Kinetic Power 

In Mil, we showed that Fkin is an important param- 
eter in classifying RL AGN. We have selected from the 
2LAC a subset of 152 blazars (as well as the detected 
radio galaxies) with estimates of the 300 MHz isotropic 
lobe emission, which is scaled to estimate Fkin as in Mil. 
For 46 sources with adequate coverage of the high-energy 
SED, we estimated F p jcfrom a two-sided parabolic fit 
to the 2LAC data in combination with X-ray data taken 
from NED 6 and/or Swift/ BAT (Cusumano et al. 2010). 
For an additional 106 sources, F p ,ic is scaled from the 
rest-frame luminosity at 1 GeV estimated from the 2LAC 
fitted spectrum: log F Pj ic = log FiGeV + v( a ~ l) 2 , where 
77 = 1.5 if a > 1, 77 = 5 otherwise. The difference between 
the two estimates is typically less than 0.2 decades in log 
L. 

We use the radio core dominance R ce = log (F core /F ext ) 
as a tracer of the orientation angle, where the core lu- 
minosity is measured at 1.4 GHz and the extended at 
300 MHz (see Mil). As discussed in Mil, the absolute 
normalization between orientation angle and R ce depends 
on Fkin, with R ce decreasing with increasing orientation 
angle for a given Fkin- 

In Figure 2 we show the 145 sources which have a 
known redshift. F Pi ic is positively correlated with both 
F ce and Fkin- The OLS bisector fits to the combined 
blazar-RG sample shown have slopes of 1.1, 1.5, 1.4, 1.3, 
1.5 (±0.1) from lowest to highest bin in Fkin- The corre- 
lation between F p j c and R ce in each group is significant 
and positive (R values from 0.37 to 0.89); the apparent 
gamma-ray output of a blazar is therefore a strong func- 
tion of both Fkin and the orientation angle. 

2.4. Emission Mechanisms and Their Beaming Patterns 

For both synchrotron and SSC, the beaming pattern is 
F = L'S p+a (Dernier 1995), where F assumes isotropic 
emission in the galaxy frame, F' is the solid-angle in- 
tegrated luminosity in the comoving jet frame, S is the 
Doppler factor, and a the spectral index. For the expo- 
nent p+ct, p=3 for a ‘moving blob’, or p=2 for a station- 
ary feature in a continuous jet (Lind & Blandford 1985), 

6 NASA Extragalactic Database: http://ned.ipac.caltech.edu/ 
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TABLE 1 

Predicted Correlation Slope Values 


index values 
Pr Psyn P 7 

SSC 

L'y Rce Rp 

- Rce 

EC 

L ' y Rce Rp Rce 

3 

3 

3 

1. 1-1.3 


0 

1.7-2 

0. 6-0.7 

2 

3 

3 

1 . 6-2 


0 

2.4-3 

0 . 8-1 

2 

2 

3 

1 . 6-2 


0 

2.4-3 

1 . 2 - 1 . 5 

2 

2 

2 

1 . 2 - 1 . 5 


0 

1.7-2. 5 

0 . 8-1 


with the possibility of different values for different emit- 
ting regions. For the EC case, the pattern is Lec = 
L' EC S p+1+2a (Dernier 1995; Georganopoulos et al. 2001). 


The slope of T Pl ic versus R ce depends on the gamma- 
ray emission process. From the synchrotron beaming, 

Rce = fo§(-^core/ A ex t) = {(Pr T Gy) log(5 + Cl, (1) 

where the factors c n depend on the unbeamed luminosity 
and jet power, but do not affect the slope. We take a = 
1 for the peaks, and eliminating 6 we obtain: 

log Tpeak — ( j ) Rce H“ C 2 , (2) 

\Pr +a r J 

where b = 1 + p 7 (SSC) or b = 3 + p 7 (EC). For ref- 
erence we list in Table 1 the expected slopes for various 
cases of p values assuming typical values of a r (0—0.5). 
The slopes in Figure 2 appear to be more consistent with 
SSC under a p 7 , p r = 3 or p 7 , p r = 2 scenario. However, 
the lack of Fermi-detected misaligned FR2 sources (par- 
ticularly, with Rce < —0.5) is likely to affect the slope 
for the two highest bins in Lkinhr Figure 2; we show for 
reference the upper limits to L p jc 7 for four RG in the 
Mil sample. We take the rightmost of these points as 
the most constraining: including it in the fit for the log 
Tkin — 45—45.5 bin increases the slope to 1.5±0.2, and 
for the highest bin to 1.7±0.3. These slopes, and our as- 
sumption that the Lorentz factor (T) of the gamma-ray 
emitting plasma is similar to that in radio, are discussed 
below. 

3. THE COMPTON DOMINANCE IN FERMI BLAZARS 
3.1. EC in powerful FSRQ 

For sources with a high-energy component dominated 
by SSC, the Compton dominance ( R v ) should remain 
constant over all orientations (no correlation between 
i? p and R ce expected) . For the EC case, 

R p = log(L E c/Ls ) = (p 7 - p syn + 2) logi) + c 3 . (3) 

Using Equation 1 we have the relation 

R p = ( Pt - Psyn + 2 \ Rce +CA . (4) 

As shown in Table 1, we expect a correlation with a slope 
from 0.6— 1.5 under the simplest assumptions for sources 
emitting gamma-rays by EC. 

Figure 3 shows R p versus R ce for the sample discussed 
in Section 2.3. Sources are divided into broad bins of 

' From the flux limit (> 100 MeV): log /ioo = —8.3 + 
(5/6 )(q— 1.2), with average a =1.3 for our sample and L p ,ic scaled 
from LiQeV (Section 2.3). 



Fig. 2. — Estimated IC peak luminosity (L Pl ic) versus radio core 
dominance (Rce) for 145 blazars and ten RG (circled) detected 
by Fermi. The correlation between L p ,ic and Rce becomes clear 
when binning by n . OLS bisector fits to sources grouped by 
Lkin have slopes 1.1, 1.5, 1.4, 1.3, 1.5 (±1; Pearson’s R values 0.81, 
0.89, 0.82, 0.37, 0.44) beginning with the lowest bin. For an EC 
emission model, these slopes are predicted to be higher (> 1.7); 
however, for the high-power sources, the lack of Fermi detections 
for misaligned (i.e., low Rce) sources will affect the apparent slope. 

Lkin- These bins include a mix of sources in terms of op- 
tical type (FSRQ or BL Lac); however all the very high 
power (VHP, Lkin > 10 45 5 ) are ‘strong-jet’ types in our 
classification due to their Lkin and position in the syn- 
chrotron plane, including many apparent BL Lacs, which 
likely suffer dilution of their broad lines by the jet emis- 
sion (Georganopoulos & Marscher 1998; Ghisellini et al. 
2011). For the moderate-power sources, the division into 
strong/weak roughly follows the FSRQ/BL Lac divide. 

No trend between R p and R ce is evident for blazars 
overall. However, when VHP sources are selected, a pos- 
itive correlation emerges. We also show the upper lim- 
its on l? p for VHP sources from the sample of Mil not 
detected by Fermi 8 . VHP sources from the 2LAC and 
Mil with upper limits on L ext 9 are also shown; because 
of the binning on Lki n which is scaled from L ext , we also 
show as connected gray points the maximum R ce such 
that Lkin > 10 45 ' 5 erg s _1 . The OLS bisector fit through 
all VHP points (including upper/lower limits) gives a 
slope of 1.1±0.1; this is a lower limit since most of the 
R v upper limits are on the lower half of the correlation. A 
sub-sample of sources with R ce values contemporaneous 
to Fermi , calculated from the average 15 GHz 10 flux over 
the time-frame of the 2LAC as measured by the Owens 
Valley Radio Observatory (OVRO, Richards et al. 2011), 
gives an identical slope. 

The non-zero slope suggests EC emission in VHP 
sources; however, when taken together the slopes in Fig- 

8 See previous footnote, with index a estimated from the relation 
for detected sources: a=— 0.2xi/ p!ayn + 3.9 

9 Estimated at 300 MHz from the lowest-frequency SED point 
with spectral index a = 1 . 2 . 

10 scaled to 1.4 GHz by —1.2 in log L 


4 


Meyer, Fossati, Georganopoulos, Lister 



R ce [1.4 GHz/300 MHz] 


Fig. 3. — The Compton dominance ( R p ) versus radio core dom- 
inance (itce)- A positive correlation is seen for the most powerful 
sources (shaded), as predicted for EC emission. The slope of the 
OLS bisector fit to the shaded points is l.liO.l, though the pres- 
ence of the upper limits on R p and lower limits on R ce (range shown 
by black lines) makes this a lower limit on the real slope. Weak 
jets (black points) as well as moderately powerful FSRQ (open 
triangles) show no trend, as expected for simple SSC. 

ures 2 and 3 are inconsistent with a simple EC scenario, 
as seen from Table 1. We have assumed that 5 is the same 
for all energies; if the radio-emitting plasma is slower, 
then the expected slopes above will increase from the 
quoted ranges for both figures. Assuming truly higher 
slopes for the VHP sources in Figure 2, concordance can 
be achieved with a general p= 3 model with a slower T in 
the radio. 

3.2. The Test of Superluminal Motions 

As shown by Lister & Marscher (1997), a substantial 
number of sources are expected to be seen at orientation 
angles smaller than that of 6 = 1/T which maximizes 
the superluminal speed j3 app ■ We therefore expect the 
highest /3 app values at R p ( R ce ) values less than the max- 
imum. Indeed if the highest R p (R ce ) sources are most 
aligned (within 1/T), their /3 app values should be rela- 
tively small (noting that these may be under-represented 
as their /3 app are difficult to measure). In Figure 4 (upper 
panel) the VHP sources are plotted as in Figure 3, col- 
ored according to /3 app . The highest /? app are observed 
in the middle of the correlation, as expected. This is 
also seen in the lower panel, where we have estimated 
the maximum /l app in bins of width 0.25 in R ce , running 
parallel to the correlation line. 11 



-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 


Rp 

Fig. 4. — Upper: The same plot as Figure 3, for high-power 
sources (Z/kin > 10 45,5 erg s — 1 ) with measurements of maximum 
apparent jet speed (/3 app ). Values are taken from Jorstad et al. 
(2001, 2005); Kellermann et al. (2004); Lister et al. (2009), and 
Kharb et al. (2010). The appearance of low /3 app sources with high 
core and Compton dominance supports the interpretation that 
these sources are aligned very close to the line of sight. Lower: 
Sources are binned along the shown correlation (running perpen- 
dicular to the dashed line in units of 0.25 of R ce & s illustrated by 
the two dotted lines), and the maximum /3 app is estimated. The 
increase in /3 app , max from lower R C e/Rp towards a maximum, fol- 
lowed by a decrease towards the highest values is qualitatively as 
expected if R ce and R p increase with alignment. 

Lkin strong jets are not. However, the typical error in log 
Lkinof 0.7 is a significant factor, as it may be the case 
that the restriction Lki n >10 45 ' 5 erg s _1 is simply that 
which is high enough to avoid any contamination with 
SSC sources. Thus the true boundary Lkin for the SSC to 
EC transition is probably lower than 10 45 5 (indeed, a few 
FSRQ of moderate Lkin are co-spatial in Figure 3 with 
the VHP sources). However, our findings imply that (a) 
VHP sources are dominated by EC, and (b) many lower 
Lkin strong jets appear at a similar (low) R p regardless of 
-R ce , suggesting that EC is not important for some part 
of the population. 

Which emission mechanism dominates is related to 
how fast the flow is: assuming the gamma-ray emission of 
strong jets is inside the BLR or MT, the comoving energy 
density of the external photon field is r 2 f 7 ext . The syn- 
chrotron energy density is U s = L s /(47rc 3 t 2 al .r 6 ), where 
t va r is the observed variability timescale in hours, L s the 
synchrotron luminosity, and we have assumed d=F. EC 
will dominate over SSC provided T 2 U ex t > U s , or for T 
greater than a transition value 


3.3. Why is EC Only Apparent in Powerful Jets? 

What is apparently depicted in Figure 3 is that 
strong jets at high Lki n are EC dominated, while lower 

11 The estimator for /3 a pp,max is derived from the order statis- 
tic >; = niax(,3 app ) for a bin of size n using the estimator for 
the unknown upper bound of a uniform distribution /3 a pp.max = 
W(n + l)/n with variance /3 ap p !max /(n 2 + 2n) (e.g. Casella 2001). 


r tr = 16.2 


L s /10 47 

(Car/6) 2 (Lext/10" 4 ) 


(5) 


Typical estimates of (Text in the BLR and MT differ by 
~ 100 (Cext.BLR ~ 2.6 x 10 -2 erg cm -3 and Lext.MT ~ 
2.6 x 10 -4 erg cm -3 ; Ghisellini et al. 2009b). In the MT 
we have Ttr.MT ~ 14.4, higher by a factor of 100 1 / 8 = 1.8 
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than in the BLR. where Ptr.BLR ~ 8.1. 

However, VLBI studies indicate that T > 10 for most 
FSRQ (e.g. Jorstad et al. 2005); if the GeV emission site 
is in the BLR, it is difficult to explain the lack of EC sig- 
nature for some strong jets. It has also been found that 
sources with higher Lkin produce on the average faster su- 
perluminal motions (Kharb et al. 2010). The connection 
between Lkin and T, and an emitting region in the MT , 
can then explain a transition to EC at high Lkin : Strong 
jets at lower Lkinare also slower and as long as T < P tr 
they are SSC emitters, exhibiting a Compton dominance 
independent of radio core dominance. As Akin increases, 
r also increases and once T > P tr , the photon held of 
the MT begins to dominate, producing a Compton dom- 
inance that increases with increasing radio core domi- 
nance, as seen in Figure 3. If we adopt longer t v ar , both 
Ttr,BLR and T^mt decrease, requiring that essentially 
all strong jets are EC emitters, which our data does not 
support. On the other hand, with a shorter f var ~1 hr 
(see e.g. Foschini et al. 2011) we obtain T tr = 12.7 for the 
BLR and T tr = 22.5 for the MT, making either location 
plausible. 

4. CONCLUSIONS 

The gamma-ray luminosities of Fermi-detected radio- 
loud AGN have been shown to depend strongly on both 
Akin and the orientation angle. We find the first collective 
evidence for external Compton emission in high-power 
jets (Akin > 10 45 ' 5 erg s^ 1 ): as can be seen in Figure 3, 
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